Effect of impurities on transport through organic self-assembled molecular films from 

first principles 
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We calculate electron transport through molecular monolayers of saturated alkanes with point 
defects from first principles. Single defects (incorporated Au ions, kinks, dangling bonds) produce 
deep localized levels in the gap between occupied and unoccupied molecular levels. Single defects 
produce steps on the I-V curve, whereas pairs of (unlike and like) defects give negative differential 
resistance peaks. The results are discussed in relation to the observed unusual transport behavior 
of organic monolayers and compared with transport through conjugated polythiophenes. 

PACS numbers: 85.65.+h 



Recently, an interest towards the development and ap- 
plication of Langmuir-Blodgett films and self-assembled 
monolayers (SAM) has risen dramatically, which is 
in part related to widespread efforts in the area of 
molecular-size devices and sensors. At the same time, 
the fundamental understanding of these films is assessed 
as "very limited" [1]. Further, it has been concluded 
some decades ago that the conduction through absorbed 
[2] and Langmuir-Blodgett [3] monolayers of fatty acids 
(CH2)„ (we shall customarily refer to them as Cn) is 
associated with defects. In particular, Polymeropoulos 
and Sagiv have studied a variety of absorbed monolayers 
from C7 to C23 on AI/AI2O3 substrates and found that 
the exponential dependence on the length of the molec- 
ular chains is only observed below the liquid nitrogen 
temperature of 77K, and no discernible length depen- 
dence was observed at higher temperatures [2]. The cur- 
rent varied strongly with the temperature for LB films 
on AI/AI2O3 substrates in He atmosphere (believed to 
hinder the AI2O3 growth) [3], which is not compatible 
with elastic tunneling. The 100- fold increase of resistance 
of the annealed films was attributed to the "annealing" 
of "defects" with time. An absence of direct tunneling 
through self-assembled monolayers of C12-C18 has been 
reported by Boulas et al. [4]. On the other hand, an 
exponential dependence of monolayer resistance on the 
chain length L, Ra- oc exp{l3a-L), [5,6] and no tempera- 
ture dependence of the conductance in C8-C16 molecules 
was observed over the temperatures T = 80 — 300K, a sig- 
nature of the elastic tunneling [6] . 

The electrons in alkane molecules are tightly bound to 
the C atoms by tr— bonds, and the band gap (between 
the highest occupied molecular orbital, HOMO, and 
lowest unoccupied molecular orbital, LUMO) is large, 
~ 9 — lOeV [4], Fig. la. In conjugated systems with 
TT— electrons the molecular orbitals are extended, and the 
HOMO-LUMO gap is correspondingly smaller, as in e.g. 
polythiophenes, where the resistance was also found to 
scale exponentially with the length of the chain, i?7r oc 



exp(/3^L), with (3^ = O.35A-1 instead of /S^ = LOSA"! 
[5]. In contrast with the temperature- independent tun- 
neling results for SAMs [6], recent extensive studies of 
electron transport through 2.8 nm thick eicosanoic acid 
(C20) LB monolayers at temperatures 2K-300K have es- 
tablished that the current is practically temperature in- 
dependent below T < 60K, but very strongly tempera- 
ture dependent at higher temperatures T — 60 — 300K 
[7]. 

A large amount of effort went into characterizing the 
organic thin films and possible defects there [1,8,9]. It 
has been found that the electrode material, like gold, 
gets into the body of the film, leading to the possibil- 
ity of metal ions existing in the film as single impurities 
and clusters. Electronic states on these impurity ions are 
available for the resonant tunneling of carriers in very 
thin films (or hopping in thicker films, a crossover be- 
tween the regimes depending on the thickness). Depend- 
ing on the density of the impurity states, with increasing 
film thickness the tunneling will be assisted by impurity 
"chains" , with an increasing number of equidistant im- 
purities [10]. One-impurity channels produce steps on 
the I-V curve but no temperature dependence, whereas 
the inelastic tunneling through pairs of impurities at low 
temperatures defines the temperature dependence of the 
film conductance, G(T) oc T'^^^, and the voltage depen- 
dence of current I{V) oc V'^^'^ [11]. This is a generic be- 
havior observed in inorganic systems with defects [12,13]. 
The tunneling may be accompanied by interaction with 
vibrons on the molecule, causing step-like features on the 
I-V curve [14,15]. During processing, especially top elec- 
trode deposition, small clusters of the electrode material 
may form in the organic film, causing Coulomb block- 
ade, which also can show up as steps on the I-V curve. 
It has long been known that a strong applied field can 
cause localized damage to thin films, presumably due to 
electromigration and the formation of conducting fila- 
ments [16]. Recent spatial mapping of the conductance 
in LB monolayers of fatic acids with the use of conducting 
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FIG. 1. Local density of states and transmission as a func- 
tion of energy for (a) C13 with Au impurity and (b) C13 with 
Au impurity and H vacancy (dangling bond) . Middle sections 
show closeups of the resonant peaks due to deep defect levels 
with respect to the HOMO and LUMO molecular states. The 
HOMO-LUMO gap is about 10 eV. 



AFM has revealed damage areas 30-lOOnm in diameter 
appearing after a "soft" electrical breakdown, sometimes 
accompanied by a strong temperatm'c dependence of the 
conductance [17]. A crossover from tunneling at low tem- 
peratures to an activation-like dependence at higher tem- 
peratures is expected for electron transport through or- 
ganic molecular films. There are recent reports about 
such a crossover in individual molecules like the 2nm long 
Tour wire with a small activation energy Ea ~ ISOmeV 
[18]. Our present results suggest that this may be a re- 
sult of interplay between the drastic renormalization of 
the electronic structure of the molecule in contact with 
electrodes, and disorder in the film, illustrated by Fig. 2 
(inset). 

We report the ab-initio calculations of point-defect as- 
sisted tunneling through alkanedithiols S(CH2),iS and 
thiophene T3 (three rings SC4) self-assembled on gold 
electrodes. The length of the alkane chain was in the 
range n = 9 — 15. We have studied single and double 
defects in the film: (i) single Au impurity. Figs, la, 2a, 
(ii) Au impurity and H vacancy (dangling bond) on the 
chain. Fig. lb, 2c, (iii) a pair of Au impurities. Fig. 2b, 
(iv) Au and a "kink" on the chain (one C=C bond instead 
of a C-C bond). Single defect states result in steps on 
the associated I-V curve, whereas molecules in the pres- 
ence of two defects generally exhibit a negative differen- 



tial resistance (NDR) . Both types of behavior are generic 
and may be relevant to some observed unusual transport 
characteristics of SAMs and LB films [2,3,7,17-19]. 

We have used an ab-initio approach that combines 
the Keldysh non-equilibrium Green's function (NEGF) 
method with self-consistent pseudopotential-based real 
space density functional theory (DFT) for systems with 
open boundary conditions provided by semi-infinite elec- 
trodes under external bias voltage [20,21]. To con- 
struct the system Hamiltonian, we use self-consistent 
DFT within the local-density approximation. The atomic 
cores are defined using standard norm- conserving nonlo- 
cal pseudopotentials [22] , and we expand the Kohn-Sham 
wave functions in a fireball [23] s-, p-, d- real-space atomic 
orbital basis [20] . We have checked that the results are re- 
liable with respect to the size of the basis set. The Greens 
function is determined by direct matrix inversion. Exter- 
nal bias is incorporated into the Hartree potential, and 
in this way the nonequilibrium charge density is iterated 
to self-consistency with imposed global charge neutrality, 
after which the current- voltage (I-V) characteristics have 
been computed. All present structures have been relaxed 
with the Gaussian98 code prior to transport calculations 
[24] . The conductance of the system at a given energy is 
found from 

g{E, V) = ^Tr [(ImE,,0 Gj^, (ImS,,,) G^j] , (1) 

where G^^^^ is the retarded (advanced) Green's function, 
S is the self-energy part connecting left (1) and right (r) 
electrodes [20], and the current is found as 



1=^1 dE[f{E~qV/2)-J{E- 



■qV/2)]g{E,V). (2) 



The equilibrium position of an Au impurity is about 3A 
away from the alkane chain, which is a typical Van-der- 
Waals distance. Note that the conductance of alkanes is 
dominated by tunneling through the wide HOMO-LUMO 
gap (see below) and, therefore, is not sensitive to fine de- 
tails of molecule-electrode coupling, unlike in the case of 
conjugated molecules, where the sensitivity to the geom- 
etry may be very high [25]. As the density maps show 
(Fig. 2), there is an appreciable hybridization between 
the s- and d-states of Au and the sp- states of the car- 
bohydride chain. Furthermore, the Au ion produces a 
Coulomb center trapping a 6s electron state at an energy 
Ci — — 0.35eV with respect to the Fermi level, almost 
in the middle of the HOMO-LUMO - lOeV gap in C„ 
(DFT gaps are usually smaller than in experiment). The 
tunneling evanescent resonant state is a superposition of 
the HOMO and LUMO molecular orbitals. Those or- 
bitals have a very complex spatial structure, reflected in 
an asymmetric line shape for the transmission. Since the 
impurity levels are very deep, they may be understood 
within the model of "short-range impurity potential" 
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FIG. 2. Current-voltage characteristics of an alkane chain 
C13 with (a) single Au impurity (6s-state), (b) two Au impu- 
rities (5d and 6s-states on left and right ions, respectively), 
and (c) Au impurity and H vacancy (dangling bond). Double 
defects produce the negative differential resistance peaks (b) 
and (c). Inset shows the density of states, transmission, and 
stick model for polythiophene T3. There is significant trans- 
mission at the Fermi level, suggesting an ohmic I-V charac- 
teristic for T3 connected to gold electrodes. Disorder in the 
film may localize states close to the Fermi level (schematically 
marked by arrow), which may assist in hole hopping transport 
with an apparently very low activation energy (0.01-0. leV), 
as is observed. 

[26]. Indeed, the impurity wave function outside of the 
narrow well can be fairly approximated as 

(p{r) — C ex.p{—Kr)/r, (3) 

where C — \/ k/2tt, k is the inverse radius of the state, 
?i^K^/2m* = Ei, where Ei = A — is the depth 
of the impurity level with respect to the LUMO, and 
A=LUMO-i^ is the distance between the LUMO and 
the Fermi level of gold and, consequently, the radius of 
the impurity state 1/k is small. The energy distance is 
A « 4.8eV in alkane chains (CH2)„ [4] (w 5eV from 
DFT calculations), and m* ^ 0.4 is the effective electron 
tunneling mass in alkanes [6] . For one impurity in a rect- 
angular tunnel barrier [26] we obtain the Breit-Wigner 
form of conductance g{E) 

^^^^^^(i.-.f + (r. + r.)V4' 

where TL(R)/'f'- is the tunneling rate from impurity to 
the left (right) electrode, Vl{R) ^ / D = Foe-^'^-^K^) , 



where t is the tunneling amplitude from impurity to the 
electrode, D is the electron band width in the electrodes. 
Using this model, we may estimate for an Au impurity in 
C13 {L = 10.9A) the width F^ = Ffl = 1.2 x 10-^ which 
is within an order of magnitude compared with the value 
1.85 X lO^^eV calculated from (1). The transmission is 
maximal and equals unity when E = and F^, = F^, 
which corresponds to a symmetrical position for the im- 
purity with respect to the electrodes. The electronic 
structure of the alkane backbone, through which the elec- 
tron tunnels to an electrode, shows up in the asymmetric 
lineshape, which is substantially non-Lorentzian, Fig. 1. 
The current remains small until the bias has aligned the 
impurity level with the Fermi level of the electrodes, re- 
sulting in a step in the current , Ii ~ ^Tqc^'^^ (Fig. 2a). 
This step can be observed only when the impurity level 
is not very far from the Fermi level F, such that biasing 
the contact can produce alignment before a breakdown 
of the device may occur. 

The most interesting situations that we have found re- 
late to the pairs of point defects in the film. Indeed, some 
impurity levels on two sites may move in and out of res- 
onance, see Eq.(6), and this obviously will lead to peaks 
on the I-V curve, or, equivalently, to Negative Differen- 
tial Resistance (NDR) . If the concentration of defects is 
c <^ 1, the relative number of configurations with pairs 
of impurities will be small, oc c^. However, they give an 
exponentially larger contribution to the current. Indeed, 
the optimal position of two impurities is symmetrical, a 
distance L/2 apart, with current I2 oc e^'^^l'^ . The con- 
ductance of a two-impurity chain is [26] 

(p. 4g2 TLTRth 

Qi2[E) = — o- 

^^\{E-ei+irL){E-e2 + tTR)-tj,f 

(5) 

For a pair of impurities with slightly differing energies 

ti2 ~ '2{Ei+E2)e If. T^jj^gj-g J. jg j-jj^g (fistance between 

them. The interpretation of the two-impurity channel 
conductance (5) is fairly straightforward: if there were 
no coupling to the electrodes, i.e. F^ = F^, = 0, the 
poles of gi2 would coincide with the bonding and an- 
tibonding levels of the two-impurity "molecule" . The 
coupling to the electrodes gives them a finite width and 
produces, generally, two peaks in conductance, whose rel- 
ative positions in energy change with the bias. The same 
consideration is valid for longer chains too, and gives an 
intuitive picture of the formation of the impurity "band" 
of states. The maximal conductivity 912 = q'^ /irh occurs 
when ei — 62, ^l^r — ti2 — where F2 is the width 
of the two-impurity resonance, and it corresponds to the 
symmetrical position of the impurities along the normal 
to the contacts separated by a distance equal to half of 
the molecule length, ri2 = L/2. Under bias voltage the 
impurity levels shift as 
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= Cio + qVzi/L, (6) 

where Zi are the positions of the impurity atoms counted 
from the center of the molecule. Due to disorder in the 
film, under bias voltage the levels will be moving in and 
out of resonance, thus producing NDR peaks on the I-V 
curve. The most pronounced negative differential resis- 
tance is presented by a gold impurity next to a Cn chain 
with an H-vacancy on one site, Fig. lb (the defect corre- 
sponds to a dangling bond). The defects result in two res- 
onant peaks in transmission. Fig. lb. Surprisingly, the H 
vacancy (dangling bond) has an energy very close to the 
electrode Fermi level F, with e, = — O.leV (Fig. lb, right 
peak in the middle). The relative positions of the reso- 
nant peaks move with an external bias and have crossed 
at 1.2V, producing a pronounced NDR peak in the I-V 
curve, Fig. 2c. No NDR peak is seen in the case of an Au 
impurity and a kink C=C on the chain bc;c;aiisc the energy 
of the kink level is far from that of the Au 6s impurity 
level. The calculated values of the peak current through 
the molecules were large: Ip « 90 n A/molecule for an Au 
impurity with H vacancy, and w 5 nA/molecule for dou- 
ble Au impurities. We have observed a new mechanism 
for the NDR peak in a situation with two Au impurities 
in the film. Namely, Au ions produce two sets of deep 
impurity levels in Cn films, one stemming from the 6s or- 
bital, another from the 5d shell, as clearly seen in Fig. 2b 
(inset). The 5d-states are separated in energy from 6s, 
so that now the tunneling through s-d pairs of states is 
allowed in addition to s-s tunneling. Since the 5d-states 
are at a lower energy than the s-state, the d- and s-states 
on different Au ions will be aligned at a certain bias. 
Due to the different angular character of those orbitals, 
the tunneling between the s-state on the first impurity 
and a d-state on another impurity will be described by 
the hopping integral analogous to the Slater-Koster sda 
integral. The peak current in that case is smaller than 
for the pair Au-H vacancy, where the overlap is of ssa 
type (cf. Figs. 2b and 2c). Thiophene molecules behave 
very differently since the tt— states there are conjugated 
and, consequently, the HOMO-LUMO gap is much nar- 
rower, just below 2 eV. The tail of the HOMO state in the 
T3 molecule (with three rings) has a significant presence 
near the electrode Fermi level, resulting in a practically 
"metallic" density of states and hence ohmic I-V charac- 
teristic. This behavior is quite robust and is in apparent 
disagreement with experiment, where the tunneling has 
been observed [5]. However, in actual thiophene devices 
the contact between the molecule and electrodes is ob- 
viously very poor, and it may lead to unusual current 
paths and temperature dependence [19]. 

We have presented the first parameter-free DFT cal- 
culations of a class of organic molecular chains incorpo- 
rating single or double point defects. The results suggest 
that the present generic defects produce deep impurity 
levels in the film and cause a resonant tunneling of elec- 



trons through the film, strongly dependent on the type of 
defects. Thus, a missing hydrogen produces a level (dan- 
gling bond) with an energy very close to the Fermi level 
of the gold electrodes F. In the case of a single impurity, 
it produces steps on the I-V curve when one electrode's 
Fermi level aligns with the impurity level under a certain 
bias voltage. 

The two-defect case is much richer, since in this case we 
generally see a formation of the negative-differential re- 
sistance peaks. We found that the Au atom together with 
the hydrogen vacancy (dangling bond) produces the most 
pronounced NDR peak at a bias of 1.2V in C13. Other 
pairs of defects do not produce such spectacular NDR 
peaks. A short range impurity potential model repro- 
duces the data very well, although the actual lineshape 
is different. 

There is a remaining question of what may cause the 
strong temperature dependence of conductance in "sim- 
ple" organic films like [CH2]„. The activation-like con- 
ductance oc exp{—Ea/T) has been reported with a small 
activation energy Eg, ~ 100 — 200meV in alkanes [7,18] 
and even smaller, 10-100 meV, in polythiophenes [19]. 
This is much smaller than the value calculated here for 
alkanes and expected from electrical and optical measure- 
ments on Cn molecules, i5a ^ A ^ 4eV [4]. which cor- 
respond nicely to the present results. In conjugated sys- 
tems, however, there may be rather natural explanation 
of small activation energies. Indeed, the HOMO in T3 
polythiophene on gold is dramatically broadened, shifted 
to higher energies and has a considerable weight at the 
Fermi level. The upward shift of the HOMO is just a con- 
sequence of the work function difference between gold and 
the molecule. In the presence of (inevitable) disorder in 
the film some of the electronic states on the molecules will 
be localized in the vicinity of Ep. Those states will assist 
the thermally activated hopping of holes within a range 
of small activation energies ^0.1 eV. Similar behavior is 
expected for Tour wires [18], where Ep — HOMO ^ 1 eV 
[20] (c), if the electrode- molecule contact is poor, as is ex- 
pected. Metallic protrusions (filaments), emerging due 
to clcctromigration in a very strong electric field, and 
metallic, hydroxyl, etc. inclusions [16,17] may result 
in a much smaller tunneling distance d for the carriers 
and the image charge lowering of the barrier. The im- 
age charge lowering of the barrier in a gap of width d 
is AU — q^\n4/{ed), meaning that a decrease of about 
3.5eV may only happen in an unrealistically narrow gap 
d = 2 - 3A in a film with dielectric constant e = 2.5, 
but it will add to the barrier lowering. More detailed 
characterization and theoretical studies along these lines 
may help to resolve this very unusual behavior. We note 
that such a mechanism cannot explain the crossover with 
temperature from tunneling to hopping reported for sin- 
gle molecular measurements, which has to be a property 
of the device, not a single molecule [18] . The work has 
been supported by DARPA. 
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